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Abstract 



We analyze the evolution of the perturbations in the infiaton field and metric following the 
end of infiation. We present accurate analytic approximations for the perturbations, showing 
that the coherent oscillations of the post-infiationary condensate necessarily break down long 
before any current phenomenological constraints require the universe to become radiation 
dominated. Further, the breakdown occurs on length-scales equivalent to the comoving post- 
infiationary horizon size. This work has implications for both the infiationary "matching" 
problem, and the possible generation of a stochastic gravitational wave background in the 
post-infiationary universe. 



1. Introduction 



A complete inflationary scenario must provide a "graceful exit" from the phase of acceler- 
ated expansion and (re)-thermalize the universe, setting the stage for a hot big bang. In 
many models this process is efficiently driven by parametric resonance sourced by the oscil- 
lating infiaton field [T|2|3] . The rich phenomenology of preheating is often contrasted with 
"standard" reheating, where the universe is rethermalized via the decay of the infiaton into 
standard matter. Given that the couplings responsible for the decay must be small in or- 
der to protect the infiaton potential against loop corrections, reheating is assumed to be a 
slow and gradual process (see e.g. ^). Without parametric resonance, it might appear that 
reheating could take place almost arbitrarily slowly: we must only insist that the universe 
reheat to MeV scale temperatures, a scale about eighteen orders of magnitude below the 
characteristic energy at the end of GUT scale infiation, in order to be rethermalized in time 
for Big Bang Nucleosynthesis and to generate a cosmological neutrino background [5]. 

In standard single field scenarios the infiaton rolls slowly through the infiationary portion 
of the potential, and then performs rapid, damped oscillations about a minimum of the 
potential after infiation is over. The frequency of these oscillations is much higher than the 
expansion rate of the universe, so p is almost constant during a single oscillation. Conversely, 
the pressure p fiuctuates rapidly: at maximum amplitude = and p = —p = V{ipmax) while 
at the bottom of the potential V{ip) = and p = p = 0^/2. The equation of state parameter, 
w, thus varies between ±1 during the oscillatory phase. Eventually, the amplitude of these 
oscillations will be small enough so that the potential is well approximated by the first 
term in its Taylor expansion around the minimum, i.e a quadratic potential. The average 
equation of state can then be shown to be zero. This leads to the conclusion that in the 
post-infiationary universe the energy density decays like and the scale factor grows like 
a{t) oc t^/^, mimicking a matter dominated universe. Since the energy density can decrease 
by 72 orders of magnitude before reheating must occur, the universe can expand by a factor 
of up to lO^'^ during this effectively matter dominated periodic] 

In this paper, we carefully analyze the evolution of the perturbations in both the scalar 
field and the metric during this period of coherent oscillations. Given that there are two 
well-separated time scales in this system (the Hubble time and the oscillation period), we 
can obtain accurate analytic approximations for the field evolution by expanding in the ratio 
of the time scales. In the limit that the infiaton is coupled only to gravity and has a purely 
quadratic potential, we show that the perturbations grow and the system becomes nonlinear 
after the scale factor a{t) has increased by a factor of about 10^. We expect that keeping 
higher order terms in the Taylor series expansion of the potential about its minimum lead to 
an even earlier onset of non-linearities, making this a conservative estimate. Consequently, 
while we are free to assume that the universe is matter dominated from the end of infiation 
to the MeV scale, the phase of coherent oscillations is necessarily much shorter. We see that 
the perturbations that become non-linear first are those for which the physical momentum 
becomes comparable to the expansion rate of the universe at the end of infiation. The 
coherent oscillations thus fragment on length scales roughly equal to the (comoving) horizon 
size at the end of infiation. 

^This is a conservative limit but the precise scale one adopts here has no impact on our conclusions. 
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We can understand our results heuristically by noting that in a universe dominated by 
pressureless dust, sub-horizon perturbations grow hnearly with a(t). Given the apparent 
red tilt of the primordial perturbation spectrum, perturbations on scales near the post- 
inflationary horizon volume have a lower initial amplitude than those at longer scales, but 
they also spend more time inside the horizon during this effective matter dominated era. 
The growth of these modes during this period more than compensates for the diminished 
initial amplitude, and shorter modes become nonlinear before longer modes. However, modes 
whose frequency is higher than the oscillation frequency of the coherent fleld "resolve" the 
oscillations. At this point our the analogy with a simple p = fluid breaks down, and modes 
with higher frequencies will be suppressed. 

Given the age and the simplicity of this model, these issues have been touched upon a 
number of times in the past. Nambu and Sasaki consider a related problem in the context of 
the invisible axion [6], while Nambu and Taruya |7] write down the equations of motion for 
the inflaton perturbations and scalar metric fluctuations, showing that they are described 
by a Mathieu-like equation. Leach and Liddle [S] compute the spectrum of perturbations for 
modes that leave the horizon near the end of inflation, while AssaduUahi and Wands [QflO] 
show that nonlinearities formed during a generic matter dominated phase likely leads to 
gravitational wave production, and a high-frequency, stochastic background of gravitational 
waves in the present-day universe. Recently, Jedamzik, Lemoine and Martin [H] (see also 
[12] ) discussed the gravitational growth of perturbations in this system, highlighting the 
parametric resonance "hidden" in the scalar fleld perturbations (see also [7]). 

The purpose of this analysis is to understand and clarify the connection between the 
detailed evolution of inflaton and metric perturbations. In particular, we carefully explore 
the basis of the simple analogy between the post-inflationary coherent oscillations and a 
dust-dominated universe. We show that the universe necessarily becomes inhomogeneous 
(and the perturbations nonlinear) if reheating is delayed long enough, and thus compute the 
maximum duration of a phase of coherent oscillations. We are careful to perform our post- 
inflationary analysis using quantities that are well-deflned at all times: some perturbation 
variables which are well-deflned during inflation have the background fleld velocity ip in the 
denominator, and thus become briefly singular during each oscillation of the fleld. The use 
of variables that remain smooth is particularly helpful for the numerical calculations we use 
to verify our analytic results. 

The layout of this paper is as follows. In Section 2, we briefly review linear perturbation 
theory for a single scalar fleld with canonical kinetic term minimally coupled to gravity. We 
focus on the appropriate choice of variables for this analysis, since the standard Mukhanov- 
Sasaki variable exhibits singularities for modes inside the horizon during the oscillatory 
phase. In Section 3, we study the evolution of perturbations that are near-horizon sized at 
the end of inflation for a quadratic potential. We calculate an appropriately deflned power 
spectrum, recovering the results by Leach and Liddle |8] where our calculation overlaps with 
theirs. Even though the underlying equation of state is rapidly oscillating, we flnd that the 
average growth of the density contrast of the modes is linear in the scale factor once they 
re-enter the horizon just like one would expect for a matter dominated universe. Finally, we 
compute when the flrst nonlinearities appear, verifying the estimate given above. We see 
that the flrst nonlinearities are found on scales comparable to the horizon size at the end 
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of inflation. In Section 4, we discuss our results and highlight a number of open questions 
identified by our analysis. 



2. General Preliminaries 

We consider single field models of inflation in which the inflaton has a canonical kinetic term 
and is minimally coupled to gravity. We first review the description and evolution of the 
perturbations in order to identify the most appropriate variables (see [13] for a thorough 
introduction to cosmological perturbation theory). In our conventions the action is 



S = [ dl^x^/^R + / d'^xy/—g 



16nG 



(2.1) 



Working to first order in the perturbations, we make the usual decomposition into a position 
independent background and inhomogeneous perturbations: 

(/9(t,x) = (^(t)+5(^(t,x). (2.2) 

Ignoring tensor modes, we fix the gauge so that the line element is 

ds^ = -N'^dt^ + a^{t)5ij{dx' + N'dt){dx^ + N^dt) . (2.3) 

The background geometry is governed by the usual Friedmann equation 



while the homogeneous scalar field obeys 

^+2Hip + V'{i^) = Q, (2.5) 

where H = d/a is the Hubble parameter and a prime indicates a derivative with respect to 
the field. 

The quadratic part of the action will govern the evolution of the perturbations. As usual, 
we use the Hamiltonian and momentum constraints to express the perturbations in the lapse 
and the shift, 6N = N — 1 and A^*, in terms of the scalar field perturbations. At first order, 
the momentum constraint implies 

6N = ^S^. (2.6) 
H 

Taking this result and the Hamiltonian constraint we see 

8^N^ = -47rG|,| (S:5^] . (2.7) 
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Using these expressions for SN and A^*, the background equations of motion, and integrating 
by parts, the quadratic part of the action for the perturbations becomes 



S = ^ I dtd^x a^{t) 



2r_2 0,2 17-2 r_2 



(2.9) 



where 

e = -— = 47rG^ and ^ = ^. 

are the Hubble slow-roll parameters. 

As long as ^ 7^ - which is certainly true during inflation - we can eliminate V"{^^ 
using the equation of motion. To be more specific, one can differentiate equation (12. 5p with 
respect to time, and obtain an expression for TpV"{Tp). Provided that ^ 7^ 0, we can divide 
by Tp^ integrate by parts, and drop the boundary term to obtain the action 



S 



dtd^x a^(t) 



6(p riHSip 



(2.10) 



Varying this action with respect to 5ip yields the Mukhanov- Sasaki equation, which can be 
brought into its standard form by making the identification 5ip = v/a, and using conformal 
time r = dt'/a{t') as the independent variable. 

Here we are interested in the transition between inflation and non-accelerated expansion, 
and the post-inflationary field oscillations. At the turning points of these oscillations, = 
and consequently e = 0, so that rj becomes singular. These singularities arise because the 
boundary term that was dropped in eliminating V"{^) becomes infinite when ^ = 0. We thus 
choose not to eliminate V"{lp) from the action, and instead express the slow- roll parameters 
in terms of the potential and its derivatives. Using the identities 



(3 - e)H^ = SnGVC^) 



and 



the action becomes 



(6 + r/ - 2e)eH^ 



H 



(2.11) 



(2.12) 



S=^ I dtd^x a^{t) 



5ip^--{d5ipf-V"{-ip)5ip' 



-167rG|\/'(<^) V - (8vrG)2|^\/(^) V 



(2.13) 



The equation of motion for the perturbations in the scalar field is then 



(5y3 + 3H5ip - ^V^Sip + V"{lp)Sip + 16nG^V'(Tp)5ip + {87rGf^V(jp)5ip = . (2.14) 



^2 



H 
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Given the invariance of this equation under spatial translations and anticipating the need to 
canonically quantize the perturbations, we look for solutions of the form 



(2.15) 



where k is the magnitude of the comoving momentum k. In our conventions the creation 
and annihilation operators a^(k') and a(k) satisfy 



[a(k),at(k')] =(27r)35(k + k'). 



(2.16) 



The modes S(fk{t) satisfy equation (12.141) with replaced by —k"^. At early times, the 
gradients dominate over the terms involving the potential. The mode equation can then be 
solved in the WKB approximation, and we can canonical quantize the field to find 



1 



exp 



-ik 



dt' 



2.17) 



a{t)y/2k 

Finally, we assume that the universe is in the Bunch-Davies vacuumjl defined by 

a(k)|0)=0. (2.18) 

Now consider modes outside the horizon, with k/aH ^ 1. In this limit the gauge 
invariant quantities TZ and ( approach constants. In our gauged] 



n, = -H 



The perturbation to the energy density of the scalar field 6pk is given by 



with 



In the limit k / aH — 



AnGTp 



H 



H 



(2.19) 

(2.20) 
(2.21) 
(2.22) 



^Departures from a Bunch-Davies initial state have been considered in the context of trans-Planckian 
corrections to the perturbation spectrum [14|15j . and would need to be substantial in order to undermine 
the conclusions reached here. 

■^Thc careful reader may be worried about the appearance of Ip in the denominator of these expressions. 
This is an artifact of the definition of these quantities at linear order in perturbations: they are only gauge 
invariant as long as ^ 7^ 0. 
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is a solution of fl2.14p for constant 7?.^°\ For this solution one has TZk = and the 

superscript (o) denotes the value outside the horizon. With one solution known, we obtain 
the second solution, 



SMt) = -C.j^J dt'[^^) (2.23) 

t. 

During an inflationary period, the second term decays exponentially in time and TZk rapidly 
approaches the constant 7?.^°\ and remains fixed until the mode re-enters the horizon. In 
the post-inflationary phase, the second solution again decays, but like rather than 
exponentially. From these two equations, we see that 

t. 

The solutions f l2.22p and fl2.23p are exact solutions of equation fl2.14p in the limit of k/aH — )■ 
and no assumption about the potential, or the existence of an inflationary phase is necessary. 
Finally, equations f l2.19p . fl2.20p . and f l2.2ip together with equation fl2.5p imply that for a 
single scalar field minimally coupled to gravity 

so that (k and TZk become equal after horizon exit when the second mode has decayed. 

In the preceding discussion, we recapitulated the usual treatment of cosmological pertur- 
bations in scalar field theories, highlighting two important points. First, we work directly 
with 6ipk which remains finite and smooth even for ^ — >■ 0, both when manipulating the 
action (by refraining from substituting the equation of motion for the V" term) and in our 
choice of the independent variable itself. The equations of motion are then both easy to 
solve numerically and convenient to analyze analytically with WKB methods. Second, we 
explicitly identify the decaying solution for 6ipk- One might expect it to be important for 
modes whose wavelength is close to the post-inflationary horizon size. Knowledge of the 
decaying solution allows us to check that it is negligible even for these modes at sufficiently 
late times. 



3. Post-Inflationary Dynamics 

We now turn to the study of the post-inflationary period, during which the fleld undergoes 
damped oscillations around the minimum of its potential. For small enough amplitude, 
almost any potential is well approximated by the leading term in its Taylor expansioi]0 

\/(^) = ^mV- (3.1) 

^Without loss of generality, we can stipulate that the minimum occurs at = 0, and we are implicitly 
setting the vacuum energy to zero. 
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Our main interest is to understand when and on what scales the first nonhnearities arise. 
Including higher order terms in the Taylor expansion turns on self-couplings of the scalar 
field and will necessarily modify the results. However, we expect that these interactions will 
generically speed up, rather than impede, the formation of non-linearities. Neglecting all 
self-couplings is thus expected to provide the most conservative estimate of the duration of 
the phase of coherent oscillations. For simplicity, we will assume that the potential is not 
only quadratic during the post-inflationary era but also during inflation. Nothing in our 
analysis depends directly on the detailed form of the perturbation spectrum at astrophysical 
scales, but we do tacitly assume that inflation ends near the GUT scale, so our analysis 
below would need minor modifications when applied to low scale inflation. 

For our numerical calculations, we assume that the pivot scale /c* = 0.002 Mpc~^ exits 
the horizon 60 e-folds before the end of inflation, and we take m = 6.35 x 10~^Mp, to be 
roughly consistent with the normalization of the power spectrum given in [5]. The precise 
value of m does not make a material difference to any of the statements below. However, 
it is worth noting that if we do have a long matter dominated phase and m need to be 
solved self-consistently. 

3.1. Background Dynamics 

By construction e = 1 when a = 0, and sets the endpoint of inflation. At this instant 



where Mp = 1/ vSttG is the reduced Planck mass. To a good approximation (p^nd ~ 



m 



Hcnd ~ ■ (3.3) 

The Hubble parameter decreases rapidly thereafter so that H <^m/2 during the oscillatory 
phase. Thus a at late times we can solve for Ip and H in the WKB approximation. At 
leading order we recover the well known result [T6] 



pit) ^ V6M„— sinfmt + A) , (3.4) 
m 

where A is some constant phase. With appropriate choice of the zero of time the expansion 
rate of the universe during this period is given by 

H^H^ = ^^. (3.5) 

The index m indicates that this is the expansion rate expected for a matter dominated 
universe. 

We use (k as a gauge invariant measure of the growth of the fractional energy density. It 
can be found by calculating TZ and using equation (12.251) . We choose the somewhat longer 
route using equations (12.191) together with (I2.20p and (I2.2ip . and verify that the two agree. 
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Figure 1: The left plot shows the number of e- folds after the end of inflation as a function of 
time. The right plot shows the Hubble constant as a function of the number of e-folds after 
the end of inflation. The orange lines represent the result of a numerical integration. The 
red curve with the longer dashes represent the analytic results of equation (13. 6p . The black 
line with the shorter dashes represents the leading order WKB approximation. 



Since the denominator in (k in equation fl2.19p is of order H^, we need to know 6pk to order 
H^. This requires us to compute the background at second order in -ffm/''^|fl 



ip{t) = V6Mp 
Hit) = H, 



H„ 



m 



3 H 

sin(mt + A) + -— cos{mt + A) 



m 



1 _ -:^sin (2(mt + A)) + 
4 m 



(3.6) 
(3.7) 



where the dots stand for corrections higher order in i?rn/''^l§ For the background scalar 
field the leading corrections have frequency m and 3m, and are suppressed by {Hm/fnY- 
The leading correction to the expansion rate has frequency 2m and is also suppressed by 
(Hm/m)'^. The leading correction to the scale factor is of order (Hm/m)'^ (see also [IE]). We 
will thus take a{t) = am{t) oc t'^^'^. These results are shown together with the leading order 
result and the result of numerical integration in Figures [1] and [2j 



3.2. Evolution of the perturbations 

We now turn to the post-inflationary evolution of the perturbations, and focus on the be- 
havior of non-relativistic modes, i.e. those with k/am ^ 1 - that is the physical wavenum- 
ber is much less than the inflaton mass. Since a{t) is monotonic, once a mode becomes 

^ This may be achieved either by going to the next order in the WKB approximation, or alternatively 
by looking for a solution schematically of the form if = 2y^2/3Mp/mt + ^ aki{mt)^^ sin((2fc — l)mt + A) + 
X;6fe/(mt)-'cos((2fc- l)mi + A) and H = 2/M + J2 Ckiimt)~' sm{{2k)mt + A) + J2 dki{mt)-' cos(2kmt + A) 
for positive integers k and / satisfying I > 2, and solving for the coefficients order by order in an expansion 
in 1/mt. This expansion is possibly asymptotic, but we work at large mt and truncate the series after a few 
terms so issues of absolute convergence do not concern us. 

^The coefficient of the correction to the field depends on terms of third order in H,n/m in our expansion 
for the background. Since we will not use these terms, we will not give them here. 
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Figure 2: These plots show the evolution of the background field ^{t) and its derivative lp{t) 
as a function of the number of e-folds after the end of infiation. The orange curves are the 
result of numerical integration. The red curve with the longer dashes represent the analytic 
results of equation (13. 6p . The black line with the shorter dashes represents the leading order 
WKB approximation. 



non-relativistic it remains so forever, in contrast to its transition across the horizon. Con- 
sequently, modes that leave the horizon during infiation are non-relativistic throughout 
the post-infiationary era. A mode which has not left the horizon as infiation ends has 
k/ttcnd > Hcnd and becomes non-relativistic after infiation is complete. From the previous 
section, the dominant contribution for these modes in the limit k/aH — )■ is given by (12. 22 p . 
Using (13. 6p . this becomes 



6Mt) = -VQMjMr^ 



cos{mt + A) 



3 H„ 



2 m 



sin(mt + A) 



3 H„ 



m 



sin(3(mt + A)) 



(3. 



where TZ^^^ is the value of TZk at times late enough so that not only Hmjra <^ 1, but also 
{k/amHm)^Hm/m <^ 1. We will see below how the second inequality arises. When the 
infiaton oscillates around its minimum, the decaying solution is suppressed relative to this 
solution by a factor of {Hm/fny, and we neglect it. This implies that all modes eventually 
approach the dominant solution and become real whether they exit the horizon during the 
infiationary period or not. 

For non-relativistic modes, we can calculate the momentum dependence to leading order 
in an expansion in k/am. For this purpose, it is sufficient to know the equation of motion 
for the perturbations (I2.14p to linear order in our expansion in H^/m. This implies that the 
last term in equation (I2.14p is negligible, and the equation takes the form 



5(pk + 3Hrn6(pk + -i^Sifk + m 



1 + 6— sin(2(mt + A)) 
m 



6(pk 



0, 



(3.9) 
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Calculating the leading correction due to k/am to leading order in Hm/m, one finds 



cos(mt + A) — sin(mt + A) + -— ^ sin(3(mt + A)) 



2 m 



1 



k 



\ drn. H„ 



m 



m 



sin(mt + A) 



. (3.10) 



Notice that the momentum dependent piece decays like t^^^ so that the correction is always 
small at late times. The amplitude of these modes thus approaches a constant, —y/QMpTl^i^\ 
This solution together with a numerical solution of equation fl2.14l) is shown for different 
values of comoving momentum in Figure [31 
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Figure 3: These plots show the evolution of the modes S(fk{t) for different values of comoving 
momentum (in inverse meters) as a function of the number of e-folds after the end of inflation. 
The orange curves are the result of numerical integration. The red curve with the longer 
dashes represent the analytic results of equation 03.81) . The black line with the shorter dashes 
represents the leading order WKB approximation. 
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This solution immediately allows us to calculate TZk and (k- One finds 



7^fc = n 



(SAl) 



The dots indicate that we have dropped terms that are of the same order in k/am but 
higher order in H^/m, as well as terms higher order in both expansions. At late times TZk 
approaches a constant, whether the modes are outside the horizon or not. Note that while 
the late time limit TZ^i^^ is well-defined, the approach to this limit is not monotonic, thanks 
to the presence of the tan(mt + A) term, ensuring that TZk acquires a "spike" whenever this 
tangent diverges. The prefactor to this term scales like t~^^^, so the spikes get narrower, but 
are always present at this order in perturbation theory. It is because of these divergences 
that we chose to work with 6ipk, which is a smooth and well-behaved function - greatly 
facilitating both our analytical and numerical discussions. 

Let us pause to understand the "microscopic" origin of the constancy of TZk- On time 
scales long compared to 1/m but short compared to 1/H, an inspection of equation (12.141) 
reveals that the modes are in parametric resonance, provided k/ttm ^ V^^Hmm, or equiva- 
lently {k / aHmY Hm/ Sm < 1, as already pointed out in [TIITT] . It turns out that this resonance 
actually yields the solution with constant TZk- To substantiate this, consider a mode that is 
observable in the CMB. At early times, the mode is inside the horizon and 6ipk{t) oscillates 
with a frequency set by its momentum. The amplitude of S(pk{t) decays like 1/a, as in equa- 
tion (12.171) . Once the mode exits the horizon (long before the end of inflation), it freezes out 
and remains constant until the expansion rate drops below the mass. At this point Sipk{t) 
becomes oscillatory, even if it is still far outside the horizon. In the absence of the term 
driving the resonance identified already in [71ITT] , the mode would decay like 1/t during this 
period, violating the constancy of 71^ outside the horizon. 

Parametric resonance is usually synonymous with exponential growth, but in this case 
the time-dependent amplitude of the term driving the parametric resonance provides just 
enough amplification to keep TZk constant for all modes with {k/ amHmYHm/'i'rrt < 1 whether 
inside or outside the horizon. Thus this "parametric resonance" ensures that not only the 
background but also the perturbations behave very similar to the ones in a matter dominated 
universe, in which IZk, or equivalently the Newtonian potential, is also constant both for 
modes inside and outside the horizon. 

Now let us turn to Cfc- From a direct calculation, we find that 



where the dots again indicate terms that are higher order in the expansions in Hm/fn-, k/am. 
This satisfies the relation (I2.25P between Ck and TZk- Further, the second term grows linearly 
in the scale factor, just as in a matter dominated universe. One can make the correspondence 
even more precise by considering a time averaged version of C,k defined as 




(3.12) 



{^Pk)T 



(3.13) 
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where the relevant quantities are averaged over some time T much longer than 1/m, but short 
compared to the expansion rate of the universe. Using equations fl2.20p . f l2.2ip and ( 13. Sp . 
one finds 



(oo) 
k 



1 



k 



15 V (Xm.Hrr 



(3.14) 



which is precisely what one expects for a matter dominated universe. Alternatively, we can 
calculate the density contrast 5k directly, where 



5 = — 



P 



_ 1^2 1 2—2 



It takes the form 



5k = 7^^) 



6 cos^(mt + A) 



5 V Clm.Hn 



(3.15) 



(3.16) 



The second term in brackets increases with time, so we expect the system to become non- 
linear, and study this in detail in 13.31 

We finish this subsection by computing For modes that exit the horizon before 

the end of infiation, 'R!'^^ = 7^^°^ In general, 7^^°^ has to be calculated numerically, but for 
modes which exit the horizon while the field is slowly rolling, the slow-roll approximation 
can be used and one has 



7^ 



H\tk 



(3.17) 



4e(4)M/ fc3 ' 

where tk denotes the time at which the mode k exits the horizon. 

We can also calculate this quantity for modes that never leave the horizon and satisfy 
k/aH ^ 1 at all times. For these modes, as long as k/am ^ y/SmH^, the solution to 
equation f l2.14p is well approximated by the WKB solution 



a3/2(t)v/2E(t) 



exp 



dt'E{t') 



with E{t) 



a{t) 



(3.18) 



Once k/ttm becomes comparable to y/3mKm, parametric resonance sets in and the WKB 
solution is no longer valid. As in [17], in this regime we look for a solution of the form 



6ipk{t) 



A{t) 



_g -j(mt+A) _^ 



Dropping terms suppressed by k/am and introducing the independent variable 

k"^ Hm 



(3.19) 



X 



(3.20) 



12 



the equation of motion for the perturbations of the scalar f l2.14p turns into a coupled system 
of equations for A and -60 

A'-iA + 3— = 0, (3.21) 

X 

B' + iB + 3- = 0, (3.22) 

X 

where the prime indicates a derivative with respect to x. This system can be solved by 
converting it into a second order differential equation for A. The initial conditions are 
chosen by matching to the WKB solution in the regime where both solutions are valid, i.e. 
for ymH <^ k/ a <^ m. One finds 

B{x) = e»Y^ - A - ^ V (3.24) 
^ ' \2x x^ 2x^ J ^ ^ 



At late times, when x <^ 1 



15 a^/'^m?H^ 

5^k -1=1 \\ cos(mt + A) . (3.25) 



From this, we obtain v}^'^ 



A;6 



9 v^A:6 ^^^^^-^^i^I^- ^^-2^^ 



Notice that H^a^ is a constant during this period, and we can evaluate it at any convenient 
time. The key result here is that the power spectrum decreases like ~ l//c^ for these modes. 

We now know the power spectrum both for modes that are highly non-relativistic near 
the end of inflation, and for modes that are highly relativistic at that time. For the transition 
region in which modes are neither highly relativistic nor highly non-relativistic at the end 
of inflation, we will rely on a numerical evaluation. The result is shown in Figure HI At 
either extreme our result recovers the slow-roll approximation for the power spectrum, and 
our (13.261) . We find excellent agreement between our analytic and numerical results in their 
appropriate regimes of validity. 



3.3. Growth of structures and the breakdown of linear perturbation theory 

For modes inside the horizon, grows linearly with the scale factor, indicating that pertur- 
bation theory will eventually break down. Let us now make this more precise. 

^Dropping terms suppressed by k/am is equivalent to treating ^ m and \B/B\ ^ in and working 

to leading order in the WKB approximation. 
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Figure 4: This plot shows the power spectrum A^(A;) as a function of comoving momentum. 
The orange curve is the result of numerical integration. The black line with the short dashes 
represents the result of the slow- roll approximation fl3.17p . The black line with longer dashes 
represents the result fl3.26p . Beyond the green vertical line modes no longer exit the horizon. 
As a consequence, one can no longer evaluate fl3.17l) . 



To ensure the validity of perturbation theory, 6N must remain small compared to the 
unperturbed value of the lapse, i.e. unity. To assess this, it is convenient to look at the 
variance 

^A^ik) . (3.27) 

The amplitude of this oscillatory function is constant in time, and thus does not signal a 
breakdown of perturbation theory at late times. The integral is sensitive to the contributions 
of modes with very long wavelengths and naively diverges for a power spectrum that is scale 
invariant in the limit k 0. However, we only expect the spectrum to be scale invariant for 
modes that exited the horizon during the inflationary period, and we expect the power in 
modes with even longer wavelengths to be suppressed due to causality, making the integral 
convergent. Under these circumstances, the integral is sensitive to the total number of e-folds 
of inflation. Since the amplitude of the observed perturbations is rather small, this does not 
pose a serious bound on the total number of e-folds, and the linear treatment is valid. 

Next, let us look at A^*. To assess its importance, we compare the contribution to the 
extrinsic curvature of the spatial slices due to A^' to the contribution from the background. 
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2.5x10"'^ - 




0.5 1 1.5 2 2.5 3 

k/m-i 



Figure 5: This plot shows the result of a numerical calculation of the integrand in equa- 
tion f l3.29p . k^A'^{k), as a function of comoving momentum. For small k the curve increases 
like k^, for large k it decays like 1/k^. The integrand peaks for values of k comparable to 
the comoving horizon size at the end of inflation. 



and require that the ratio remain small. This translates into 

1 



^^{{m\t,0))')^l. (3.28) 
Using (12. 7p . together with the results in the last subsection, one finds 

i,((a,;V<(,0)n = l(-^)7*PAiW. (3.29) 

The integrand is shown in Figure O As one would expect, the integral receives its largest 
contribution from the modes whose physical momentum at the end of inflation is close to 
the Hubble constant at that time. These are modes in the transition region where the power 
spectrum changes from nearly scale invariant to decaying like 1/k^. We see that this quantity 
grows like during the post- inflationary era, and this term is the one that fixes the onset 
of nonlinearity. The integral has excellent convergence properties both for small and large 
momenta, and can easily be performed numerically. We finclfl 

{{m\t, 0))2) ^ 1.4 X 10-12 e^^ , (3.30) 



^This formula is only valid after the relevant modes have approached the solution with TZk- This happens 
about 4 e-folds after the end of inflation. 
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where denotes the number of e-folds after the end of inflation. We conclude that the onset 
of non-hnearity occurs roughly 13 e-folds after the end of inflation. Mapping between H{t) 
and a(t), we see that non-linearities develop when H ^ 3 x 10~* -f^cnd or H ^ 10~^^Mp. If 
the universe thermalized instantaneously at this point, we would have a reheat temperature 

4 X IQii GeV 

Trh ^ —J-, , (3.31) 

where A/'eg is the effective number of relativistic degrees of freedom. 

It may be interesting to point out that variance of the density contrast 6{t, x) = 6p{t, x) /p 
approaches unity around the same time. At late times it is given by 



{{6{t, 0))2) = 36 cos^(mt + A) j dk A^(A;) 

cos2(mt + A) ( — \—] I dk kA'^{k) 



5 



4 . . .4 



^Y,[^J y^^^^A^(^). (3.32) 

The first term is small for any realistic total number of e-folds of inflation. The third term 
becomes large around 13 e-folds after the end of inflation, while the second term grows more 
slowly. 

Finally, let us consider the scalar perturbations. Two things have to be satisfied to justify 
treating the scalar field perturbations to linear order. First, their variance must remain sub- 
Planckian. At times late enough for the modes to have become real and to have approached 
the solution that oscillates with constant amplitude, this variance is given by 

/dk 
yA^(fc) . (3.33) 

The amplitude of this oscillatory function is again independent of time, and remains small 
at late times. 

Second, the scalar field perturbations should be small compared to the scalar field back- 
ground to ensure that the energy density stored in the background dominates over that in 
the perturbations so that the perturbations do not affect the background geometry. This 
requires {Sip(t, 0)^) -C ^^(t), which is valid as long as 

^|^A^(fc)^2.9xlO-V^«l, (3.34) 

where A^ denotes the number of e-folds after the end of inflation. This becomes of order unity 
around six e-folds after the end of inflation and suggests a breakdown of our perturbation 
theory even earlier than we had found by looking at the metric perturbations. As was 
argued in [17], however, this effect could be resummed so that the evolution can be studied 
anal5d;ically even after this happens and the system truly becomes non-linear only once the 
density contrast approaches unityl^ 



'We are very grateful to Alexei Starobinsky for pointing this out to us. 
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4. Discussion 



We have given a full and careful account of the evolution of the universe during a phase of 
"coherent oscillations" following inflation. We have worked in variables that are manifestly 
finite, and given reliable analytic approximations for the individual modes' evolution. As 
is well known, the background solution is pressure-free when averaged on time scales long 
compared to the oscillation time of the fields. All modes grow linearly with the scale factor, 
matching our expectations for a dust-dominated universe with vanishing pressure. Modes 
whose natural frequency is much lower than the frequency of the scalar field oscillations 
during the end of inflation possess a nearly scale invariant power spectrum. For modes 
whose frequency during the end of inflation exceeds that of the background field oscillations, 
the power spectrum decreases rapidly with increasing frequency. 

The breakdown of the analogy between the coherent oscillations of the infiaton and 
pressureless dust reflects the presence of two intrinsic time scales in our dynamical system, 
namely the Hubble time and the oscillation period for the scalar field condensate. Conversely, 
a dust-dominated universe has only the Hubble time. Consequently, the perturbations of the 
scalar field match those of dust only when we can "integrate out" the more rapid scalar field 
oscillations. Interestingly, the "long modes" for which the pressureless dust approximation is 
reliable can be shown to undergo parametric resonance while they are outside the horizon jTI] , 
and this resonance is the microphysical origin of the constancy of TZk, as we show via an 
explicit calculation. 

Phenomenologically, this distinction between "long" and "short" modes picks out a 
unique scale in the post-inflationary universe. At the end of inflation, modes which are 
approximately horizon-sized have a frequency similar to the oscillations of the coherent in- 
fiaton background. Solutions for modes well inside the horizon at the end of inflation decay 
in amplitude until they are redshifted to lower frequencies by the expansion of the universe. 
Conversely, modes that left the horizon well before the end of inflation are frozen until they 
re-enter the horizon. Consequently, modes that are roughly horizon sized at the end of in- 
flation undergo significantly more growth than other modes and will be the first for which 
nonlinear effects become important, and the breakdown of the coherent oscillations takes 
place on length scales fixed by the post- inflationary horizon size. 

Physically, the breakdown of coherence in the oscillating field implies that while we 
can assume a very long period of matter domination following inflation, the universe can 
expand by at most a factor ~ e^^ or so during the coherent oscillation phase. Beyond this 
point, Sp/p ~ 1, and the field is no longer homogeneous. Once this breakdown occurs, 
we need to move beyond first order perturbation theory, raising the prospect that even the 
simplest post-inflationary universe can have highly nontrivial phenomenology. This could 
include gravitational wave production Pll0fl2] . or even the formation of primordial black 
holes (e.g. [18]) if the overdensities grow without limit. In the latter case, the primordial 
black holes will be very small, since their mass will be set by the mass-energy inside a 
region whose size is equal to the comoving length scale at the end of inflation. These black 
holes would thus decay quickly via Hawking radiation, reheating the universe (e.g. [19]) via 
"non-perturbative" gravitational effects. Moreover, in addition to any gravitational waves 
generated via the nonlinear evolution of the field, gravitons Hawking radiated as the black 
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holes decay would generate a background of very high-frequency gravitational waves |20] . 

More generally, this analysis reinforces the importance of carefully exploring the physical 
processes which govern the universe between inflationary and nucleosynthesis scales. Models 
with parametric resonance are well known for their rich phenomenology and possible grav- 
itational wave background (e.g. [2T]|22II23II24|I25] . but it is becoming increasingly clear that 
even the simplest potentials can have nontrivial post-inflationary dynamics. We also wish to 
understand the overall thermal history of the post-inflationary epoch in order to accurately 
match present-day astrophysical scales to their counterparts during the inflationary era |26] . 
Consequently, while this model is almost as old as inflation itself |27], its phenomenology 
contains a number of important open questions. The analysis here provides a thorough and 
complete understanding of the post-inflationary oscillations, and we intend to pursue the 
issues raised here in future work. 
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